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ABSTRACT 


An  investigation  has  been  made  of  the  effect  of  the  isolated 
blade  row  in  a  single-rotor,  axial  flow  compressor  upon  a  circumferentially 
distorted  inlet  velocity  profile.  Flow  rate  and  chord/pitch  ratio  were 
varied  to  determine  the  influence  of  these  parameters  on  the  velocity  de- 
fect. Experimental  results  were  compared  with  two  theories,  both  based  on 
a  two-dimensional,  incompressible,  inviscid  actuator  disc  analysis.  It 
was  found  that  both  theories  underestimate  the  experimentally  observed  at- 
tenuation due  to  the  rotor,  however,  the  more  recent  analysis  by  Seidel 
gives  better  agreement  than  that  of  Rannie  and  Marble.  Seidel 's  theory 
correctly  predicts  more  attenuation  with  an  increase  in  chord/pitch  ratio. 
The  gain  in  attenuation  due  to  reduced  flow  rate  was  not  accurately  pre- 
dicted by  either  theory,  but  again  Seidel' s  prediction  is  closer  to  ex- 
periment . 
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CIRCUMFERENTIAL  DISTORTION  OF  THE  INLET  FLOW  IN  AN  AXIAL  COMPRESSOR 


1.   INTRODUCTION 

Since  modern,  high-performance  aircraft  must  operate  over  a  wide 
range  of  flight  conditions,  the  gas- turbine  power  plants  for  such  aircraft 
must  evolve  from  designs  which  attempt  to  provide  satisfactory  engine  per- 
formance  over  a  comparable  range  of  inlet  velocity  distributions.  Inlet 
velocity  distortions  in  axial  flow  power  plants  may  be  produced  by  separa- 
tion at  the  inlet  lip  when  operating  at  high  angles  of  attack  or  off- de- 
sign conditions,  or  by  separation  from  the  walls  when  the  flow  is  turned 
or  diffused  too  rapidly  within  the  inlet  duct.  Non- uniformity  of  the  in- 
let velocity  presentation  may  be  expected  to  result  in  reduction  in  com- 
pressor stall  margin,  and  non-uniform,  thus  inefficient  compression  of  the 
working  fluid. 

While  jet  engines  differ  in  their  acceptance  of  inlet  flow  dis- 
tortions, both  theory  and  investigations  of  representative  compressor  in- 
let stages  indicate  that  velocity  distortion  will  be  reduced  across  a  rotor 
because  that  portion  of  a  distorted  flow  having  an  energy  level  lower  than 
some  mean  value  is  energized  more  than  portions  having  an  energy  levejL  higher 
than  the  mean.  Further  investigation  of  the  problem  of  inlet  flow  distor- 
tion in  axial  turbo-machines  is  justified  in  order  to  formulate  design  cri- 
teria for  compressors  which  will  accept  more  severe  distortions,  and  for 
ducting  which  will  reduce  distortions  before  reaching  the  compressor. 
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Circumferential  distortions  of  the  inlet  flow,  which  appear  to 
cause  the  most  significant  reduction  in  compressor  stall  margin,  have  been 
the  subject  of  only  a  few  analytical  publications.  Notable  among  these  are 
the  work  of  Ehrich  (Reference  l),  and  that  of  Rannie  and  Marble  (Reference 
2).  In  both  of  these  investigations,  a  quasi- steady  actuator  disc  analysis 
was  applied  to  an  incompressible,  inviscid  flow  in  an  axial  turbomachine 
with  hub/tip  ratio  close  to  unity,  allowing  the  problem  to  be  treated  two- 
dimensionally . 

Seidel  (Reference  3),   in  a  more  recent  investigation,  has  com- 
pared  the  theories  of  Ehrich  and  Rannie-Marble  and  concluded  that  they  are 
the  same  for  the  case  of  an  isolated  blade  row.  He  suggests  an  improve- 
ment to  the  previous  theory  in  which  a  more  general,  unsteady  actuator 
disc  analysis  is  used.  The  assumption  from  previous  theory  of  constant 
leaving  angle  relative  to  the  rotor  is  relaxed  in  favor  of  the  Kutta  con- 
dition, and  the  introduction  of  unsteady  blade  forces  permits  investiga- 
tion of  the  effect  of  variations  in  chord/pitch  ratio. 

The  experimental  investigation  herein  described  served  to  com- 
pare the  two  existing  theories  with  the  performance  of  a  single-rotor, 
axial  flow  compressor.  A  single  inlet  screen  configuration  was  used  for 
all  tests,  and  data  were  obtained  for  four  flow  rates  at  each  of  two  values 
of  chord/pitch  ratio.  In  all  cases,  measurements  were  taken  at  four  axial 
locations; 

i)   (-  00  )   -  far  upstream  of  the  rotor 
ii)   (  1  )   -  immediately  upstream  of  the  rotor 
iii)   (  2  )   -  immediately  downstream  of  the  rotor 
iv)   (+  00  )   -  far  downstream  of  the  rotor 
Comparison  with  theory  is  presented  for  the  maximum  and  minimum  flow  rates 
for  each  of  the  two  chord/pitch  ratios. 
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2.   APPARATUS  AND  INSTRUMENTATION 

2.1  Compressor 

Figure  1  is  a  section  view  of  the  compressor  used  in  these  tests. 
The  isolated  rotor  and  a  portion  of  the  instrumentation  arrangement  are 
shown  in  Figure  2.  Essential  rotor  dimensions  were: 

Outside  diameter  2J.25  inches 

Hub/tip  ratio  O.75 

Blade  chord*  I.50  inches 

Blade  aspect  ratio  2.Qk 

Blade  tip  clearance         O.O35  inces 
The  blade  sections  were  NACA  65-(l2)l0  airfoils  with  a  stagger 
at  the  mean  radius  of  52.7  degrees  measured  from  the  axial  direction.  The 
blades  were  twisted  linearly  9*70  degrees  from  root  to  tip.  A  rotor  speed 
of  1000  rpm  was  used  in  the  tests  for  which  results  are  presented.  Inlet 
air  was  introduced  radially  through  damping  screens. 

2.2  Distortion  Screens 

The  choice  of  an  inlet  distortion  screen  arrangement  was  dictated 
by  the  desire  to  reduce  the  solidity  at  the  edges1 of  a  screen  segment,  thus 
avoiding  inlet  distortion  profiles  complicated  by  the  formation  of  concen- 
trated vortices  in  these  regions.  When  such  vortices  do  occur,  a  velocity 
distribution  similar  to  that  in  Figure  3  would  be  anticipated  from  poten- 
tial theory,  and  such  velocity  distributions  are  reported  in  the  experi- 
mental results  of  Reference  k. 
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The  most  desirable  inlet  screen  arrangement  provides  a  velocity 
profile  which  is  readily  defined  mathematically  and/ or  typical  of  distor- 
tion in  actual  machinery,  e.g.  a  square-wave  or  sinusoidal,  inlet  disturbance. » 
While  relatively  simple  to  represent  mathematically,  the  sine  wave  is  dif- 
ficult to  obtain  experimentally,  leading  to  the  choice  of  a  square  wave  pro- 
file generator.  For  generation  of  a  square -wave  profile,  the  ideal  solidity 
distribution  in  the  circumferential  direction  is  shown  in  Figure  4,  however 
such  screening  is  not  readily  available. 

Several  combinations  of  square-mesh  screens  of  different  solidity 
were  tested  to  obtain  a  profile  which  approximated  a  square  wave,  despite 
the  appearance  of  moderate,  vortex- induced  peaks.  The  final  arrangement 
consisted  pf  a  45  degree  screen  of  mesh  4x4,  wire  size  23  (A.W.G.),  on 
which  was  superimposed  symmetrically  a  36  degree  screen  of  mesh  6x6 
made  from  number  25  wire.  Three  of  the  larger  screens  provided  the  rig- 
idity for  an  annular  frame  which  is  shown  in  Figure  5*  In  order  to  achieve 
relative  motion  between  screen  and  probes,  the  system  of  screens  and  frame 
was  positioned  stepwise  through  a  complete  cycle  of  120  degrees  by  means 
of  a  reel  and  cord  arrangement  (Figure  6)  which  permitted  readings  every 
1.42  degrees  of  arc.  The  use  of  three  distortion  screen  segments  and  the 
location  of  the  probes  far  upstream  and  far  downstream  of  the  rotor  were 
dictated  by  considerations  outlined  in  Appendix  A. 

2.3  Probes 

Shown  in  Figure  7  is  a  three-hole  Fecheimer  tube  of  the  type  de- 
scribed in  Reference  5  which  was  used  to  measure  yaw  angle,  static  pres- 
sure and  total  pressure.  This  probe  has  three  0.012  inch  holes  drilled 
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into  a  0.125  inch  stock  tube  at  a  distance  of  0.25  inches  from  a  hemi- 
fspherical  tip.  The  drillings,  which  are  spaced  at  k2   degrees  of  arc,  con- 
nect with  three  hypodermic  tubes  contained  inside  the  stock  tube.  The  cen- 
ter hole  points  into  the  flow  when  the  pressures  at  the  two  outer  drillings 
are  equalized. 

The  yawing  of  the  probes  and  the  measurement  of  pressure  at  the 
four  stations  shown  in  Figure  8  were  accomplished  through  the  use  of  a 
pressure  switching  circuit  and  a  pair  of  strain  gage  pressure  transducers 
connected  to  a  D-C  calibrating  device.  This  transducer  equipment,  which 
incorporates  a  transistorized  D-C  amplifier,  was  manufactured  by  the 
Dynamic  Instrument  Company  of  Cambridge,  Massachusetts. 

Yaw  angles  were  measured  directly  with  protractors  mounted  on 
the  compressor  shell.  Initial  alignment  of  the  probes  involved  removal 
of  the  distortion  screens  to  assure  axial  flow  at  the  upstream  stations. 
Then  each  probe  in  turn  was  placed  in  the  farthest  upstream  station,  yawed 
until  the  center  hole  was  aligned  with  the  flow,  and  the  indicator  arm  set 
tj»iread  90  degrees  on  the  protractor. 

Volume  flow  of  air  through  the  compressor  was  controlled  by  the 
position  of  a  damper  valve  shown  in  Figure  1,  and  the  rpm  of  the  rotor.  A 
counter  mounted  on  the  compressor  control  panel  measured  damper  setting, 
while  an  electrical  strobotac  indicated  the  rotor  speed. 
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3-   EXPERIMENTAL  PROCEDURE 


The  highest  flow  rate  for  each  value  of  chord/pitch  ratio  was 
fixed  by  the  maximum  damper  opening  and  the  minimum  flow  rate  determined 
"by  the  onset  of  rotating  stall.  In  determining  the  lower  limit  of  damper 
setting  it  was  observed  that  rotating  stall  occurred  over  a  range  of  dam- 
per settings  depending  on  such  factors  as  local  disturbances  at  the  inlet 
bell-mouth  and.  random  noise.  As  a  consequence,  early  attempts  to  operate 
too  close  to  this  critical  range  resulted  in  runs  which  were  only  partial- 
ly  completed  when  rotating  stall  occurred.  A  low  flow  setting  was  finally 
chosen  which  permitted  continuous  readings,  and  two  intermediate  damper 
settings  between  the  high  and  low  flow  rates  were  selected  for  investiga- 
tion. Significant  parameters  for  the  maximum  and  minimum  flow  rates  at 
each  value  of  chord/pitch  ratio  are  listed  in  Table  3-1. 

TABLE  3-1 
Flow  Parameters 


Chord/Pitch 

1.05 

0.525 

U/  n  r 

•  567 

.460 

.514 

•391 

U-oo      (ft/ sec) 

53-6 

41.2 

50.6 

36.8 

Ux         (ft/ sec) 

49-5 

40.6 

46.2 

36.O 

U2         (ft/ sec) 

51.7 

43.4 

50.5 

39.9 

U+,co      (ft/sec) 

5%f+ 

^1*3 

46.1 

35.O 

v      (■• )  ■ ' 

9-3- 

15*3 

8.8 

l6vl 

wQ         (ft/ sec) 

91.2 

8&"6 

90.6 

81.8 

^b 

40,500 

32,100 

36,500 

27,400 
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For  the  maximum  value  of  chord/pitch  ratio  (2b/S  =  1.05),  kk 
blades  were  mounted  on  the  rotor  which  was  operated  at  1,000  rpm  through- 
out. The  resulting  low  air  velocity  (U  =  50  ft/ sec)  permitted  the  flow 
to  be  treated  as  incompressible. 

To  reduce  the  chord/pitch  ratio  by  one-half, 22  blades  were 
removed  alternately,  and  the  resulting  holes  in  the  rotor  rim  filled  by 
aluminum  hexagonal- head  bolts.  Again,  for  the  maximum  and  minimum  flow 
rates,  the  rotor  was  operated  at  1,000  rpm,  but,  in  order  to  duplicate 
the  intermediate  flow  rates  of  the  higher  solidity  case,  it  was  necessary 
to  adjust  the  damper  settings  slightly  and  operate  at  1,050  rpm.  An  in- 
crease in  the  stall-free  operating  range  due  to  removal  of  half  the  blades 
permitted  low  flow  rate  measurements  at  a  reduced  damper  setting. 

For  each  flow  rate,  measurements  of  flow  angle,  static  and  stagna- 
tion pressure  were  recorded  approximately  every  7*1°  degrees  of  screen  ro- 
tation except  in  regions  near  the  edges  of  the  distortion  screen.  Here 
readings  every  1.42  degrees  were  required  to  define  the  profile  more  ac- 
curately. The  decision  to  rotate  the  screen  rather  than  the  probes  was 
based  on: 

i)  the  knowledge  from  previous  experiments  with  the  same 
compressor  that  the  flow  without  distortion  screens 
was  axi symmetric; 
ii)  the  relative  simplicity  of  the  screen  rotation  mech- 
anism required  (Figure  6)  in  comparison  with  the  ex- 
tensive modification  necessary  in  order  to  rotate  the 
the  four  probes  through  a  complete  cycle  of  120  degrees. 
Measurements  were  taken  at  the  mean  geometric  radius  of  the  flow  annulus 
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to  avoid  boundary  layer  effects  and  thus  achieve  a  maximum  degree  of  two- 
dimensionality  . 

Prior  to  each  run,  ambient  pressure  and  temperature  were  re- 
corded from  a  standard  mercury  barometer  and  Fahrenheit  thermometer  re- 
spectively. Also  recorded  were  the  zero  (atmospheric)  readings  of  the 
two  pressure  transducers  in  terms  of  a  resistance  bridge  balancing  current. 
The  compressor  was  then  brought  up  to  speed  and  the  damper  set.  For  each 
data  point,  each  probe  in  turn  was  yawed  until  equalized  pressures  at  the 
outer  taps  produced  a  null  (balanced)  reading  on  the  indicating  dial,  and 
the  flow  angle  was  recorded.  The  stagnation  pressure  (center  hole)  and 
the  indicated  pressure  p^  at  one  of  the  outer  holes  were  then  measured 
by  comparison  with  atmospheric  pressure.  In  subsequent  calculations  these 

readings  were  related  to  the  free  stream  static  pressure  by  means  of  an 

P  -  Pi 
experimentally  determined  probe  calibration  constant  k  =  p  _  p  ;  the 

values  of  pQ  and  p  were  then  used  to  compute  velocity  from  V  =1—  |pQ  -  p. 


P 
Density  was  calculated  from  the  perfect  gas  relation  making  use  of  ambient 

pressure  and  temperature.  No  correction  was  made  for  stagnation  pressure 
gradient  since  previous  experimental  measurements  indicated  that  this  ef- 
fect was  negligible  .    .,«,*! 
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k.     ANALYSIS  OF  RESULTS 

k.l    Method  of  Presentation 


It  was  convenient  to  plot  the  experimental  results  in  terms 
of  the  quantities 

H  =  p/p  +  Uu  +  Vv 
B=  Uv  -  Vu 

(P-p/p 

for  ready  comparison  with  theery.  U  and  V  were  determined  from  graphi- 
cal integration  of  plots  of  local  velocity  components,  and  the  perturba- 
tion velocities  u  and  v  from  the  difference  of  the  mean  values  and  the 
local  components.  Perturbation  pressures  (p)  were  calculated  in  a  simi- 
lar manner. 

2       2       2 
Experimentally  determined  values  of  Hg/U  9    Q  ±/\J  ,    S  g/U  s 

(p  ]/u  and  (P2/u  are  presented  coincidentally  with  curves  of  the  same 

quantities  obtained  from  the  theories  of  Seidel  and  Rannie-Marble  in  Figures 

9-1  through  9-h.     No  presentation  of  results  is  given  for  the  intermediate 

flow  rates  since  the  effect  of  flow  rate  is  most  clearly  demonstrated  by 

comparison  of  the  highest  and  lowest  values  obtainable. 

k.2     Comparison  with  Theory 

A  brief  development  of  the  two  methods  of  analysis  used  for  the 
theoretical  results  is  given  in  Appendix  B  where  the  notation  of  Rannie  ; 
and  Marble  has  been  modified  in  the  case  of  ©  and  (P  to  agree  with  that 

of  Seidel.  Input  data  for  both  theories  was  obtained  from  51  equi- spaced 

2 
points  on  the  curves  of  H-j/U  shown  in  Figures  9« 
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Although  both  theories  underestimated  the  attenuation  of  the  ve- 
locity distortion  due  to  the  rotor,  it  is  apparent  that  the  Seidel  theory- 
gives  consistently  better  agreement  for  the  downstream  predicted  quantities 
R>>  O   2  and(P2«  The  Rannie  and  Marble  analysis  gives  a  more  accurate 
prediction  of  the  upstream  values     ®  ,  and  (?■,.     Agreement  in  either 
case  may  be  limited  by  the  assumption  of  inviscid  flow  and  by  the  dif- 
ficulty involved  in  locating  -0  accurately.  No  real  flow  is  entirely  with- 
out viscosity,  and,  since  the  Cauchy-Riemann  functions  ©  and  (P  vary  as 

fiX 

exp  (-  —  ),  an  error  of  1.5  inches  (the  distance  from  the  blade  row  mid- 
plane  to  the  nearest  probes,  Figure  8)  in  the  location  of  the  upstream 
or  downstream  £ero  point  means  a  reduction  in  these  quantities  to  - 
exp  (-  ^  gg  )=  .64  of  their  value  at  the  actuator  disc. 

The  quantity  H2/U  was  chosen  as  a  suitable  parameter  for  deter- 
mining the  degree  of  attenuation  of  the  velocity  distortion  since  this 
quantity  is  the  perturbation  velocity  within  a  constant.  That  is,H2/lT  = 
H+qq  /  U2  =  "+00  (1  +  V22/U2).  The  effect  of  flow  rate  on  attenuation  is 
shown  in  Figure  10,  and  the  effect  due  to  solidity  in  Figure  11.  The 
Seidel  theory  is  seen  to  correctly  predict  increased  attenuation  for  a 
higher  value  of  chord/pitch  ratio.  Since  the  Rannie-Marble  analysis  as- 
sumes constant  chord/pitch  ratio,  no  comparison  is  possible  for  that  theory. 
While  neither  theory  accurately  predicts  the  experimentally  observed  gain 
in  attenuation  with  a  reduction  in  flow  rate,  Seidel' s  theory  shows  im- 
proved agreement  with  experiment  over  that  of  Rannie  and  Marble. 

It  was  concluded  that,  in  general,  the  Seidel  theory  gives  a 
better  prediction  of  the  performance  of  a  single-rotor  axial  compressor 
with  circumferentially  distorted  inlet  velocity  profile. 
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4. 5  Applicability  fco  Theory 

To  test  the  applicability  of  the  experimental  results  to  the 
theories  considered,  three  assumptions  common  to  both  analyses  were  ex- 
amined. The  continuity  condition  assumes  U.^,  =  U,  =  Up  =  U+0o  •  ex- 
perimental results  (Table  3-l)  reveal  a  consistent  variation  of  6  to  9  per 
cent  between  the  highest  and  lowest  value,  with  high  readings  occurring 
at  stations  -  oo  and  2,  low  readings  at  1  and  +  oo  .  This  variation  may 
be  the  result  of  a  change  in  radial  distribution  of  velocity  with  axial 
location,  but  is  probably  due  to  the  influence  of  turbulence  in  the  wake 
of  the  screen  and  rotor.  Appendix  C  is  a  discussion  of  the  influence  of 
turbulence  on  pressure  measurements.  The  influence  of  turbulence  was 
presumed  least  at  +  oo  ,  and  this  value  of  U  was  used  as  the  representa- 
tive velocity  for  each  flow  rate. 

Theory  assumes  that  H(-  oo  ,  y)  =  H]_(0  -  y  y)  and  H2(0+  y  y)  = 
H(+oo^y);  Figure  12  indicates  that  H+o/U^  decays  before  reaching  the 
downstream  measuring  station,  probably  due  to  viscous  effects.  The  up- 
stream curves  agreed  more  closely,  but  no  conclusion  could  be  drawn  from 
their  comparison. 

Finally,  the  linearized  forms  of  the  equations  of  motion  and  con- 

U-oo 

tinuity  are  predicated  on  small  perturbations.  If  the  quantity  — is 

used  as  a  measure  of  perturbation  intensity,  the  third  term  in  the  ex- 

2  2 

pansion  (l+u/u)  =  l+2u/U  +  (u/u)  is  neglected  in  linear  theory.  Since 

U  -oo 
the  maximum  value  of  — - observed  experimentally  was  .35*  the  highest 

(.35)2 
value  of  the  neglected  term  in  comparison  with  the  term  retained  is   ,   . 

.175. 

It  was  concluded  that  experimental  conditions  sufficiently 
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approximated  the  three  assumed  conditions  to  make  the  tests  valid.  The 
degree  of  validity  would  he  imporved  with  correction  factors  based  on 
additional  knowledge  of  the  magnitude  of  effects  due  to  turbulence  and 
viscosity. 

5 .  CONCLUSION 

An  experimental  investigation  of  a  single- rot or  axial  flow  com- 
pressor has  determined  the  effect  of  the  isolated  blade  row  on  a  forced 
circumferential  distortion  in  the  inlet  velocity  profile.  For  the  dis- 
tortion-screen configuration  tested,  it  was  shown  that  the  velocity  de- 
fect was  attenuated  by  passage  through  the  moving  blade  row,  and  that  the 
degree  of  attenuation  increased  with  reduced  flow  rate  or  increased  chord/ 
pitch  ratio. 

Experimental  results  were  compared  with  two  theories,  both  based 
on  a  two-dimensional, incompressible,  inviscid  actuator  disc  analysis.  It 
was  found  that  both  theories  underestimate  the  attenuation  due  to  the  rotor, 
although  the  theory  of  Seidel  gives  better  agreement  with  experiment  than 
that  of  Rannie  and  Marble.  Seidel 's  theory  correctly  predicts  an  increase 
in  attenuation  for  increased  chord/pitch  ratio.  Neither  theory  accurately 
predicts  the  experimentally  observed  gain  in  attenuation  with  reduced  flow 
rate,  however  agreement  is  improved  with  the  Seidel  theory.  It  was  gen- 
erally concluded  that  the  Seidel  theory  gives  a  more  accurate  prediction 
of  isolated  rotor  performance  with  regard  to  a  circumferentially  distorted 
inlet  velocity  distribution. 

Comparison  was  made  between  experimental  results  and  three  limit- 
ing assumptions  common  to  both  analyses,  namely,  continuity,  constant  H(x,y) 
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upstream  or  downstream  of  the  rotor,  and  small  perturbations.  It  was 
concluded  that  the  experiments  were  applicable  to  the  theories  considered. 
Further  investigation  of  distorted  inlet  flows  in  axial  compres- 
sors might  provide  information  on: 

i)  radial  velocity  distribution  in  a  circumferentially 
distorted  flow  field; 
ii)  the  effect  of  multiple  blade  rows  including  stators, 
rotors,  and  inlet  guide  vanes; 
iii)  the  magnitude  of  viscous  and  turbulence  effects; 
iv)  an  easily- applied  relationship  between  attenuation 
and  mean  flow  (for  design  purposes). 
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APPENDIX  A 
Criterion  for  Inlet  Distortion  Wavelength  and  Axial  Location  of  Probes 


In  the  analysis  of  Reference  1  it  was  concluded  that  any  cir- 

2jtr 
cumferential  distortion  of  width  -rr-  will  induce  local  disturbances  at 

the  rotor.  However,  these  disturbances  will  diminish  exponentially  up- 

,-n|x| . 
stream  and  downstream  of  the  rotor,  the  decay  factor  being  exp  [ — % — ;. 

Since  three  sections  of  distortion  screen  (Figure  5)  were  used 

in  these  tests,  the  distortion  width  was  =££-  =  6.67it  inches  and  ^  =  3.3I+ 

inches.  For  the  probes  located  22.0  inches  upstream  and  downstream  from 

the  plane  of  the  rotor,  the  rotor  influence  was  reduced  by 

**B   rfSgjj)  =  0.001. 
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APPENDIX  B 


Methods  of  Analysis 


In  the  following  discussion,  that  portion  of  Reference  2  deal- 
ing with  a  circumferential  inlet  distortion  approaching  an  isolated  rotor 
has  been  extracted  (with  notational  modification)  to  indicate  the  assump- 
tions of  the  Rannie  and  Marble  solution  and  the  procedure  for  calculation 
of  the  corresponding  theoretical  curves  of  Figure  9* 

No  exposition  is  given  of  the  Seidel  solution,  however,  the 
differences  in  assumptions  are  emphasized,  and  the  sequence  of  develop- 
ment outlined  briefly. 

Bl  Theory  of  Rannie  and  Marble 

In  the  description  of  the  disturbance  flow  field,  let  x  and  y 
be  rectangular  coordinates  with  the  axis  of  y  parallel  to  the  plane  of 
the  cascade,  Figure  13.  The  average  velocity  components  and  pressure  are 
U,VandP,  all  constant  upstream  and  downstream  of  an  isolated  blade  row. 
The  disturbance  velocity  components  and  pressure  are  u,  v  and  p,  all  with 
average  values  zero.  The  disturbance  is  assumed  periodic  in  the  y-direc- 
tion,  and  of  wave  length  sufficiently  large,  compared  with  the  blade  gap, 
that  the  distortions  produced  by  individual  blade's  on  the  disturbance  pat- 
tern are  ignored. 

For  ffimii  disturbances  the  equations  of  motion  and  of  continuity 
in  linearized  form  are: 

JJ  3jL      +    V  Ist  =  -  4-  |£  (la) 
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Completely  equivalent,  but  more  illuminating  are  the  following  three 
equations: 
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The  first  of  these  is  the  condition  that  the  perturbation  of  total  pres- 
sure is  constant  along  the  mean  stream  lines.  The  second  two  equations 
show  that  the  combination [Uv  -  VuJ  (proportional  to  flow  angle  perturba- 
tion) and  the  static  pressure  divided  by  density  are  potential  functions 
and  satisfy  the  Catchy- Riemann  conditions.  For  convenience  in  the  analy- 
sis, the  following  definitions  are  introduced: 


H  - 

e  - 
(?  - 


f-    4-Ua  +  Vv 


Uv  -v 
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(5) 
(6) 
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The  following  relations  then  hold, 

H(x,y)  =  H(o,y-^x)  (8) 

and  two  functions  are  sufficient  to  define  the  disturbance  flow  com- 
pletely. 

For  a  semi- infinite  flow  field  upstream  or  downstream  of  a  blade 
row,  the  functions  can  be  given  in  simple  terms.   O  and  (P   approach  zero 
far  from  the  cascade  and  if  the  blade  row  is  at  x  =  0,  then  for  @(0,y) 
given,  (P(P,y)  can  be  found  from 


(P(o+,lj)  =  ^  Q  io +.  7})  cot  IT (u-^dli  =  6  (o+.y) 

r"         l  *,   '      \  (10) 

fP(o-.u)  -  -\<5)(o-  yi)  co"G-rr(u-7))d7)  =  @  (o-.u) 


Reference  6  contains  a  developemnt  of  this  relationship  between  comple- 
mentary functions. 

Conversely,  Q  (0+,  y)  =  -  (p   *(0+,  y)  and  0  (0-,  y)  =  (p  *(0-,  y) 
where  the  wave  length  is  unity.  ©  and  (P  may  be  found  for  all  values  of 
x  and  y  if  ©  (0,  y)  or  (p  (0,  y)  is  given.  Hence  the  disturbance  flow 
is  described  completely  by  H(0,  y)  and  @  (0,  y)  or  (P  (0,  y). 

If  the  flow  is  assumed  locally  quasi- steady,  three  conditions 
relating  to  conservation  of  mass,  momentum  and  energy  can  be  prescribed 
to  give  the  local  outlet  flow  in  terms  of  the  local  inlet  flow.  Let  U^, 
Vq_,  P-j_  and  Ug,  Vg,  P2  *>e  "the  average  velocity  components  and  pressures 
upstream  and  downstream  of  a  rotor.  Let  u-j_,  v-j_,  p^_  and  Ug,  Vg,  Pg  be  the 
perturbation  values  at  (0-,  y)  and  (0+,  y)  respectively.  Then  the  matching 
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conditions  at  the  blade  row  are  those  listed  in  Table  B-I  where,  for 
the  equations  in  H,  ©  and  (p  ,   the  following  values  are  specified  in 
any  experiment: 


it 

2 


R   =  -   ^'YWY-Y)    -|-    w-v,   +  t\j_  p  _  p 


Y_=Y  ^   m 


UY 


u 


"vp  ■•■ 


A 


UY 


y   =_  U'^  Y'-(vJfY) 
flc  UY 


The  symbol  P  represents  a  blade  loss  coefficient  which  is  a  prescribed 
function  of  tanp^.  Pr*  is  the  slope  of  the  Pr  (tanf^)  curve  at  a  point. 

TABLE  B-I 
Boundary  Conditions  for  Rannie-Marble 


Matching  Condition 


Continuity 


Constant  Leaving 
Angle 


Equation 


+  Li. 


Vr-^v,)  =  tan  Q 
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Equation  in  H,  G,  (P 
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Blade  Loss  Speci- 
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The  first  two  equations  of  Table  B-I  may  be  written  as: 

0,  -  e£  +-  Ar  S\   +■  B,  ft   -  o  <«■> 

e,  +  <*r  f;  +  /3r  (^  =  h,  w 

The  complementary  equations  to  these  two  from  the  operation  of  Eqn.  10 
are: 

<P,  +<P2  -4.G,  -  B,GS,  =0  a?) 


<P,   -  <*r ©J  +  pr St,     =   H"  W 


* 
I 


Except  for  a  phase  shift  (important  only  if  a  specific  streamline  is  de- 
sired), H(-  co   ,  y)  =  H-^Q-,  y)  is  a  given  function.  Also  known  are  the 
constants  Ay,  Br,  ctr,  £r,  and  tfr.  Equations  11  through  14  determine  (j±, 

(j^i    ®  2.   an<i  ®  2*  The  ^2  can  be  foun<i  from  the  third  matching  condi- 
tion, 

•  Hz  =  ©z    +  *r  <£  (15) 

and  the  solution  is  complete. 
B2  Theory  of  Seidel 

For  a  flow  field  prescribed  as  for  the  theory  of  Rannie  and 
Marble,  Seidel  has  chosen  to  relax  the  assumption  of  constant  leaving 
angle.  Experimental  (evidence  of  the  invalidity  of  this  assumption  is 
shown  in  Figure  Ik. 

Making  use  of  the  Kutta  condition,  Seidel  determines  the  in- 
cremental forces  on  a  single  blade  in  the  flow  field  due  to  perturba- 
tion velocity  components  normal  and  parallel  to  the  mean  relative 
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velocity.  In  this  calculation,  the  Rannie-Marble  analysis  is  used  as  a 
first  approximation  to  obtain  perturbation  components  u-^,  u^,   v-j_,  v2  which 
are  then  resolved  into  components  parallel  and  perpendicular  to  wQ.  The 
net  y- component  of  "blade  force  per  chord"  (non-dimensional)  per  unit 
span  is  then  calculated, and  this  quantity  multiplied  by  chord/pitch  ratio 
to  obtain  the  "force  per  pitch"  per  unit  span.  The  influence  of  blade  row 
effectiveness  is  introduced  as  a  lattice  coefficient  which  is  a  function 
of  chord/pitch  ratio,  stagger  and  camber. 

The  three  boundary  conditions  for  the  Seidel  analysis  are  listed 
in  Table  B-II  where  the  blade  force  specification  has  replaced  the  constant 
leaving  angle  assumption  of  Rannie  and  Marble.  Reference  3  includes  a  de- 
tailed treatment  of  the  equivalence  of  the  expressions  in  u,  v,  p  and  those 
in  H,  0  ,  (P  .  The  constant  5    values  are  fixed  by  the  mean  flow,  and 
the  values  of  H,  Q     ,  and  (?   are  computed  in  a  manner  analogous  to  that 
of  Rannie  and  Marble. 

TABLE  B-II 
Boundary  Conditions  for  Seidel 


Match  Condition 

Equation 

Equation  in  H,  ©>(P 

Continuity 

U,  =  U, 

Relative  Bernoulli 
Equation 

(n  *n  (it    2      d&   i   £ 
(id      ' 

Blade  Force  Speci- 
fication 

-aFy    zU(v,-v2)  +  u(Y-X) 

oS 

•z,  Q  -©  -  -3  (p  - 

(zo     \          *    fa     1 

^5     <-<>   i 
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Bj  Calculations 

The  calculation  time  required  for  each  theoretical  curve  of 
Figure  9  can  be  appreciated  when  it  is  remembered  that  a  system  of  four 
simultaneous  equations  must  he  solved  for  each  of  51  input  points  (Hq/u). 
To  reduce  the  manual  computation  time  for  both  methods  of  analysis,  it 
was  convenient  to  program  the  numerical  work  for  the  IBM  'JOk-  Electronic 
Data  Processing  Machine  at  the  M.I.T.  Computation  Center.  Flow  diagrams 
and  further  details  of  the  programming  and  computation  procedure  are  found 
in  Reference  5* 
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APPENDIX  C 
Influence  of  Turbulence  on  Pressure  Measurements 

In  the  measurement  of  flow  immediately  behind  the  distortion 
screen  or  the  rotor,  it  is  probable  that  pressure  readings  are  in  er- 
ror due  to  the  high  turbulence  level  in  these  regions.  The  pressure 
tap  in  turbulent  flow  no  longer  measures  only  the  pressure  of  the 
stagnant  air  just  inside  the  opening,  but  includes  a  non-zero,  time- 
average  increment  of  pressure*  resulting  from  the  turbulence  velocity 
component  normal  to  the  plane  of  the  opening. 

"Reference  7  indicates  that  the  turbulence- induced  pressure 
components  for  static  and  stagnation  measurements  are  identical  for 
isotropic  turbulence.  Thus,  for  velocity  calculations  with  isotropic 
turbulence,  the  indicated  pressures  may  be  used  without  correction, 
since  the  turbulence  components  are  self- cancelling. 

There  is,  however,  experimental  evidence  that  the  turbulence 
immediately  behind  a  grid  is  anisotropic.  In  Figure  5  and  Figure  6  of 
Reference  8,  for  a  flow  velocity  comparable  to  that  of  the  present  ex- 
periments, the  longitudinal  component  of  turbulence  is  shown  to  be  pre- 
dominant, even  at  some  distance  from  the  turbulence  producing  grid. 

In  the  experiments  for  this  thesis,  pressure  measurements  in 
the  region  immediately  downstream  of  a  distortion  screen  demonstrate  a 
measure  of  scatter  (not  observed  in  the  region  between  screens)  which 
may  be  attributed  to  turbulence.  Based  on  the  evidence  of  Reference  8, 
this  turbulence  may  be  assumed  anisotropic,  in  which  case  the  three-hole 
probe  measures  a  turbulence- induced  increment  of  stagnation  pressure  which 


is  larger  than  the  turbulence- induced  component  of  static  pressure. 

Thus  in  mean- velocity  calculations  for  low  speed  flow  (that 


is,  where  the  flow  is  assumed  incompressible  and  V  ~  / pQ  -  p  )  in  regions 
where  anisotropic  turbulence  is  indicated  (e.g.,  immediately  downstream  of 
a  distortion  screen  or  rotor)  a  value  higher  than  the  true  time-mean  value 
would  be  anticipated. 
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FI6.3       VELOCITY     DISTRIBUTION 
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FIG. 4  IDEAL      SOLIDITY      DISTRIBUTION      FOR 

A      SQUARE      WAVE      PROFILE 
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FIG.  7     THREE- HOLE  FECHEIMER  TUBE  AND 

DETAILS 
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FIG.     9-1       PERTURBATIONS    UPSTREAM    AND    DOWNSTREAM,  THEORY 
AND    EXPERIMENT.     2  b/ S  =  105,  U/Jlr  =  .562 
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FIG.     9-1      CONTD. 
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FIG.    9-3  2b/S  =525,   U/ilr=  .514 
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FIG.    9-3     CONTD. 
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FIG.  9-4  2b/S  =  .525,  U/ar  =  .39l 
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FIG.   9-4    CONTD. 
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FIG. 10    EFFECT      OF      FLOW      RATE      ON       ATTENUATION,     THEORY 
AND       EXPERIMENT.     2b/S  ■  1.05 
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FI6.II      EFFECT   OF    SOLIDITY   ON    ATTENUATION,   THEORY    AND 
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FIG.  12         COMPARISON     OF     EXPERIMENTAL      H    /U*    AND      H    /U* 
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